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ABSTRACT 

In  the  mean  field  approximation  of  nonlinear  relativistic  o-CO-  p  model,  we  study  the  influence  of  the 
density-dependent  coupling  constants,  between  nucleons  and  mesons  on  the  symmetry  energy  S(pB)  of  infinite  nuclear 
matter,  in  four  different  density-dependent  formalism.  We  find  greater  Tp  leads  to  greater  Ksym  and  L  whenT aa),  c  and 
d  are  fixed  and  indicate  larger  S(pB )  in  high  density  region.  In  addition,  the  density  dependence  of  Ya  co  p  make  S(pB ) 
smaller  in  high  density  region,  and  they  make  Ksym  and  L  more  sensitive  to  the  changing  of  parameters  at  different 
density. 
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1.  INTRODUCTION 

The  nuclear  symmetry  energy,  which  is  the  difference  of  energy  per  nucleon  between  pure  neutron  matter  and 
symmetric  nuclear  matter,  characterizes  the  isospin  dependent  part  of  the  equation  of  state  (EOS),  of  asymmetric  nuclear 
matter.  The  nuclear  symmetry  energy  is  closely  related  with  the  neutron  skin  of  a  nucleus  and  pressure  of  crust  matter  of 
neutron  star,  so  it  is  an  important  topic  in  nuclear  physics  and  astrophysics11’41.  The  information  of  the  symmetry  energy  is 
especially  interesting  in  related  fields,  and  many  theoretical  and  experimental  efforts  have  been  performed  [5"71. 
The  3H/3He  ratio  is  sensitive  to  the  nuclear  symmetry  energy  at  sub-saturation  densities,  within  the  newly  updated  version 
of  the  Ultra-relativistic  quantum  molecular  dynamics  (UrQMD)  model181.  The  impact  of  the  rearrangement  term  and 
momentum  dependence  of  the  single-particle  potential,  on  the  density  dependence  of  the  nuclear  symmetry  energy  and 
nucleon  effective  mass  has  been  studied  in  details  [91.  Both  the  isospin-singlet  and  isospin-triplet  components  of  the 
potential  energy,  play  an  important  role,  in  determining  the  symmetry  energy,  when  the  Fock  diagram  is  introduced  within 
the  covariant  density  functional  (CDF)  theory1101.  The  findings  reveal  that,  the  ratio  of  relative  yield  of  light  charged 
particles  poses  better  candidate,  to  probe  the  density  dependence  of  nuclear  symmetry  energy,  by  using  the  yield  of  various 
fragments,  in  central  collisions  of  various  isotopic  and  isobaric  colliding  pairs  [I1].  The  breakup  reactions  induced  by  the 
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polarized  deuteron  beam  at  about  100  Mevlu  ,  provide  a  more  stringent  constraint  to  the  symmetry  energy  at  sub  saturation 
densities,  via  the  calculations  on  a  novel  reorientation  effect  of  deuteron  attributed  to  is  o  vector  interaction,  in  the  nuclear 
field  of  heavy  target  nuclei [12].  An  analytic  expression  for  the  symmetry  energy  as  a  function  of  its  slope  parameter  L  is 
found,  by  using  a  modified  quark-meson  coupling  model  [13].  A  covariance  analysis  reveals  that,  the  additional  fitting 
protocol  reduces  the  uncertainties  in  the  nuclear  symmetry  energy  coefficient,  its  slope  parameter  as  well  as  the  neutron- 
skin  thickness  in  20&Pb  nucleus  by  50%[141.  Assessments  based  on  the  sensitivity  matrix  reveal  that,  the  properties  of 
extreme  neutron-rich  systems  play  a  predominant  role,  in  narrowing  down  the  uncertainties  in  the  various  symmetry  energy 
parameters  [15].  Four  quantities  deducible  from  nuclear  structure  experiments  have  been  claimed  to  correlate  to  the  slope 
parameter  L  of  the  symmetry  energy  [1  ].  Results  from  classical  molecular  dynamics  simulations  of  infinite  nuclear 
systems,  show  an  excellent  agreement  with  the  experimental  data  and  corroborate  the  claim  that,  the  formation  of  clusters 
has  a  strong  influence  on  the  symmetry  energy,  in  the  liquid-gas  coexistence  region1171.  The  slope  of  the  symmetry  energy 
plays  an  important  role,  in  determining  the  boundary  and  properties  of  the  mixed  phase,  in  the  study  of  the  liquid-gas  phase 
transition  of  stellar  matter,  with  the  inclusion  of  the  finite-size  effect  from  surface  and  Coulomb  energies1181.  Some  mean 
field  models  associated  with  similar  values  of  L  at  saturation  density,  and  pertaining  to  different  families,  yield  a  greater- 
than-expected  spread,  in  the  neutron-skin  thickness  of  the  208  Pb  nucleus1191.  And  the  main  uncertainty  of  symmetry  energy 

in  finite  nuclei  is  found  to  be  related  to  the  1 4  term,  by  investigating  the  model  dependence  of  7 4  term  in  theoretical 
symmetry  energy  for  a  few  popular  mass  models1201. 

The  recommended  value  of  symmetry  energy  at  saturation  density  is  3 1  MeV ,  and  the  behavior  of  symmetry 
energy  below  saturation  density  is  now  much  better  known,  while  the  density  above  saturation  is  still  not  yet  well 
constrained,  and  the  predictions  from  different  models  strongly  diverge.  There  are  two  different  forms  of  the  density 
dependence  of  the  symmetry  energy;  the  stiff  dependence  is  the  symmetry  energy  increases  monotonically  as  the  density 
increases,  while  the  soft  is  the  symmetry  increases  initially  up  to  normal  nuclear  density  and  then  decrease  in  higher 
density.  It  is  essential  to  constrain  the  form  of  the  density  dependence  of  the  symmetry  energy  for  understanding  the 
nucleon-nucleon  interaction  and  the  structure  of  the  compact  stellar  objects  such  as  the  neutron  star.  The  main  reason  of  the 
uncertain  symmetry  energy  is  related  to  nuclear  force  and  its  spin  and  isospin  dependence  which  have  not  been  understood 
thoroughly. 

The  specific  values  of  L  and  Ksym  are  model  dependent,  the  ref.  [21]  gives  L  =  46.7  ±  1 2. 8  MeV  and 
Ksym  =  -166.9 ±  168.3 MeV ,  the  studies  on  neutron  skin1221  gain  L  =  59 ±  1 3.0MeV  ,  astrophysical  observations  of  neutron 
star  radii  and  masses1231  recommend  L  =  48  ±4.0 MeV ,  while  our  conclusion  on  Ksym  and  L  will  be  listed  in  Section  3  in 
detail. 

In  this  paper,  we  study  the  influence  of  density-dependent  coupling  constants  on  symmetry  energy  of  nuclear 
matter  in  the  relativistic  mean  field  (RMF)  approximation  of  the  relativistic  a  -a>-  p  model.  In  Section  2  the  mean  field 
approximation  of  nonlinear  relativistic  o-co-  p  model  and  three  density-dependent  formulas  of  meson-nucleon  coupling 
will  be  introduced. 

In  Section  3  the  numerical  results  will  be  discussed  in  detail,  and  the  paper  closes  in  Section  4  with  a  summary 
and  concluding  remarks. 
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2.  BASIC  THEORIES 

The  Lagrange  density  in  nonlinear  relativistic  o-co—p  model  is 
L  =  Yj¥b a- TnY^af  -|r 'pypc-V1  )yB 

B  1 

--a^eeT  +-mlcolia>t1  +-(dModf‘eT-m2cr2)  +  Yjipx(ir/Jd,‘ -m^ 

-—pay  ■  Ppv  +  -mlp  ■  pP  ccr3  — - da 4 

Where,  i/sb(B  =  n,  p) ,  y/A ,  <x ,  p  and  co  are  the  field  operator  of  baryon  (including  neutron  and  proton)  and 
lepton,  <7  ,  co  and  p  meson,  respectively.  F^ ,  Vp  and  F„,  are  the  density-dependent  coupling  constants,  and  mn  ,  mp  ,  ma , 
mp  and  mm  are  the  mass  of  neutron,  proton,  o ,  co  and  p  meson,  respectively. 

From  this  we  can  obtain  the  energy  density  and  pressure  of  nuclear  matter. 
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The  binding  energy  of  a  nucleon  and  the  incompressibility  coefficient  are  given  by 


E(pB,S)~ - mB 

Pb 


(4) 


2  d2E{pB,S  =  0) 

0  Pbo  d2pB  pB=p‘° 


(5) 


Where,  pB  is  number  density  of  baryons,  and  pB0  is  the  saturation  density. 

The  symmetry  energy,  its  slope  and  curvature  at  saturation  density,  are  described  by 


S(PbY- 


e-P'-P" 


2  dpB 


18 


Pb  o 


(6) 

(7) 

(8) 

The  density-dependent  formulas  of  a  are  defined  as  following.  Case  I  is  introduced  from  references  |241, 
while  Case  II  and  Case  III  are  introduced  by  our  work  in  order  to  reduce  the  quantity  of  coupling  constants; 


L~3p  3S(^)| 
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K  =()n2 
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Case  I: 

r  ,(pb)  =  8,Mx)  (9> 

Mx)  =  a1  +  b,(x  +  d,)l  ,i  =  cr,  a  (10) 

‘  1  +  Cj(x  +  d ,)2 

(.Pb)  =  bp  exp[  —ap  (x  -  1)]  (H) 

where  x  =  pB  /  PBq  ,  the  eight  constants  in  (10)  are  restricted  by  /.(i)  =  i,  /"(i)=/"(i)  and  /."(0)  =  0. 

Case  II: 

rk(pB)  =  ak+iogpb‘ ,  k=a,  a,  p  (12) 

Case  III: 

rj(pB)  =  Jaj+bje(-',Bh)  ,  j  =  a,  CO ,  p  (13) 

3.  NUMERICAL  RESULTS 

In  formulas  (10)— (13),  the  constants  a br ,  ci  and  d,  =<J ,  CO  and  p),  together  with  c  and  d  in  the 

self-interaction  of  o  meson  are  adjusted  to  fit  the  binding  energy  of  nuclear  matter  E(pB0)  =  —\5.15MeV  and 
S{PBq)  =  31  MeV  at  pB0  =  0.17  jm~3 .  At  the  same  time,  the  incompressibility  coefficient  of  the  EOS  should  be  between 
220  and  300 MeV. 


The  symmetric  energy  S(pB )  in  Case  I,  Case  II  and  Case  III,  which  correspond  to  the  parameters  R(T0j\  , 

Raa>5  and  R(TOj6  in  table  1,  are  shown  in  Figure  1.  The  values  of  and  S{pB)  at  pB  =0.5 frn~3  and  the  values  of 

S(PB),  Ksym  anc*  L  at  saturation  density  pB  =  0. 1 7/w  3  are  shown  in  Table  1,  while  the  specific  values  of  other 

parameters  are  shown  in  Table  2,  several  other  models  are  present  in  Table  2  for  comparison.  As  we  can  see,  greater  values 
of  lead  to  stiffer  S(pB)  for  R<7(0\  (Case  I)  in  high  density  region  and  greater  values  of  Ksym  and  L  at  saturation 

density.  Similar  phenomena  happen  in  other  two  different  density-dependent  formulas  of  Raco5  (Case  II)  and  RaO) 6 
(Case  III). 


Impact  Factor  (JCC):  4.8764 


NAAS  Rating  3.73 


Influence  of  Density-Dependent  Coupling  Constants  on  Symmetry  Energy  of  Nuclear  Matter 


11 


Table  1:  The  Data  Corresponding  to  Figure  1 


aP 

bP 

cp 

r>(/>s=o.5) 

S(pB  =  0.5  )(MeV) 

S(pB  =0.17  )(MeV) 

UMeV) 

K^JMeV) 

RoaA-a 

0.04 

6.90 

/ 

6.39 

77.25 

31 

82.03 

-34.41 

RoaA 

RoaA-b 

0.24 

6.91 

/ 

4.33 

59.42 

31 

66.30 

-102.53 

RooA-c 

0.74 

6.91 

/ 

1.64 

46.38 

31 

26.90 

-107.22 

Ratal 

0.74 

6.94 

/ 

1.65 

47.82 

31 

26.76 

-98.24 

RotS 

0.74 

7.26 

/ 

1.73 

41.43 

31 

18.64 

-121.43 

Read 

0.74 

6.22 

/ 

1.48 

55.06 

31 

42.47 

-76.53 

Roof -a 

5.47 

0.24 

/ 

5.96 

74.33 

31 

77.85 

-36.43 

Roof 

R  ctg)5  -  b 

3.84 

0.54 

/ 

4.96 

65.57 

31 

67.75 

-57.29 

Roof  -c 

0.83 

0.74 

/ 

3.13 

53.56 

31 

49.15 

-74.31 

Roofs -a 

-64.7 

100 

-21.2 

6.14 

87.86 

31 

93.48 

19.56 

Roto 6 

Rotafs-b 

9.5 

29 

2 

5.66 

83.30 

31 

90.95 

4.79 

Roofs -c 

11.5 

50.0 

0.24 

5.01 

71.67 

31 

78.40 

-40.30 

Rcrcol 

7.56 

/ 

/ 

7.56 

80.09 

31 

81.3 

-26.4 

Table  2:  The  Coupling  Constants  of  Various  Models  (Include  Our  Work  and  Others) 


K0(MeV) 

go 

Sa> 

ao 

K 

C<T 

da 

ao 

K 

d. 

c 

d 

Racd 

259.97 

8.54 

9.85 

1.25 

0.27 

0.40 

0.91 

1.25 

0.26 

0.39 

0.92 

100 

100 

Roto! 

261.10 

8.43 

9.69 

1.25 

0.95 

1.29 

0.51 

1.24 

1.23 

1.65 

0.45 

100 

100 

Ratal 

261.60 

7.39 

8.06 

1.27 

0.26 

0.40 

0.91 

1.20 

0.60 

0.80 

0.65 

100 

100 

RaoA 

261.08 

9.68 

11.8 

1.30 

0.25 

0.40 

0.91 

1.32 

0.22 

0.37 

0.95 

0 

0 

Rato5 

260.14 

/ 

/ 

7.15 

0.39 

/ 

/ 

8.03 

0.50 

/ 

/ 

100 

100 

RacaS 

260.30 

/ 

/ 

82.8 

52.7 

0.16 

/ 

111.4 

90 

0.21 

/ 

20 

5 

Ratal 

261.64 

7.11 

6.48 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

1450 

1446 

nr  -  99 

240.00 

10.7 

13.3 

1.37 

0.23 

0.41 

0.90 

1.40 

0.17 

0.34 

0.98 

0 

0 

DD-MEl 

244.50 

10.4 

12.9 

1.39 

0.98 

1.53 

0.47 

1.39 

0.85 

1.36 

0.50 

0 

0 

HA 

233.00 

9.51 

10.4 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

7093.9 

-255.9 

NL-Bl 

280.00 

8.76 

11.8 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

15.029 

-100.9 

Figure  1:  The  Influence  of  Tp  on  the  Relation  between  the  Nuclear  Symmetry  Energy  5, Pb ,  and  pB  when  l„m ,  c 

and  d  are  Fixed.  The  Solid,  Dashed  and  Asterisk  Curves  Represent  Ratal- a,  Ratal- b  and  Ratal- c  in  Figure 
la,  Rata5-a,  Rato5-b  and  Rata5-c  in  Figure  lb,  Raco6-a,  Ratab-b  and  Raco6-c  in  Figure  lc  , 

Respectively 


In  Figure  2a,  the  value  of  S(pB)  at  saturation  density  being  fixed  on  31  MeV ,  seven  curves  with  K0  ~  260 MeV 


are  shown.  The  different  curves  represent  the  counterpart  in  Figure  2b,  Figure  2c  and  Figure  2d.  As  we  can  see, 
the  symmetry  energy  in  Case  I  is  smaller  than  the  one  in  other  three  cases  in  high  density  region  when  K0  keeps  the  same 


value  in  all  cases,  and  that  in  non-linear  model  is  the  largest.  So  the  density  dependence  of  coupling  constants  can  decrease 
the  symmetry  energy  in  high  density  region. 
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Pg/frrT 


PB/fm 

Figure  2:  The  S(pB)  of  Nuclear  Matter  with  K0  ~  260 MeV  and  Different  Values  of  Other  Parameters  (Figure  2a), 
and  the  Dependence  of  Va(op  on  Nucleon  Number  Density  (Figure  2b,  Figure  2c  and  Figure  2d).  The  Curves  in 
Figure  2a  Represent  the  Counterparts  in  Figure  2b,  Figure  2c  and  Figure  2d.  The  Curves  RaoA  —  c  and  Racol  —  4 
are  Obtained  from  Case  I  in  Table  1  and  2,  while  R<jco5-c  ,  R<JC06  -  c  and  Roofl  are  Gained  from  Case  II,  Case  III 

and  Non-Linear  Model,  Respectively 

In  Figure  3,  we  depict  the  correlations  between  L  and  K  at  saturation  density  with  K0  restricted  between  220 

and  300  MeV ,  the  distribution  range  of  L  and  Ksym  are  different  in  different  density-dependent  formulas  at  the  same 

density.  An  approximate  linear  correlation  ( Ksym  =  2.63L- 176.21)  exists  in  the  range  of  i 0 <  L<  60MeV  and 

-140  <  Ksym  <  -40 MeV  for  Case  I,  but  it  disperses  when  L  >  60 MeV ,  while  70  <  L  <  90 MeV  and  -80  <  Ksym  <  0 MeV 

for  Case  II,  40  <  L  <  HOMeV  and  -100  <  Ksym  <  -60MeV  for  Case  III,  80  <  L<  90  MeV  and  -30  <  Ksym  <  QMeV  for 

non-linear  model.  The  values  of  L  in  all  models  are  shown  in  Figure  4.  As  we  can  see,  the  range  of  the  values  of  L 
obtained  in  this  work  is  consistent  with  those  obtained  from  other  analyses. 


L  /MeV 

Figure  3:  the  Correlations  between  L  and  Ksym  Calculated  by  Different  Parameters  at  Saturation  Density  with  K0 

Restricted  between  220  and  300  MeV .  The  Asterisks,  Squares,  Triangles  and  Circles  Represent  the  Values  of  L  and 
K  in  Case  I,  Case  II,  Case  III  and  Non-Linear  Model,  Respectively 
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Case  III 

Ref.[33] 

Ref. [31]  _ 


Nonlinear 


•  Ref.[32] 


Ref.[30] 


Ref.[28] 


Ref.[29] 


Ref.[27] 

Ref.[26] 


Ref.[25] 


0  10  20  30  40  50  60  70  80  90  100  110 

L/MeV 

Figure  4:  The  Values  of  L  from  All  Models 


Similar  to  the  symmetry  energy  as  a  function  of  nucleon  density  in  Figure  1 ,  we  will  study  the  density  dependence 
of  correlation  between  L  and  Ksym  .  In  Figure  5,  the  correlation  between  L  and  Ksym  in  various  cases  at  pB  =  0.5  p0 ,  p0 , 

2p0  and  3 p0  are  shown,  respectively.  In  Case  I,  the  scatter  plot  which  is  obtained  by  changing  parameters  moves  to  the 
left  when  the  density  increases.  At  the  same  time,  the  range  of  L  is  up  to  its  maximum  around  the  saturation  density  while 
the  one  of  Ksym  keeps  minimum  as  pB  >  p0  .  On  the  other  hand,  L  and  Ksym  are  positively  correlated  when  pB  is  about 

less  than  p{) ,  and  they  become  to  be  negatively  correlated  when  pB  is  greater  than  p{] .  In  the  last  three  cases, 
similar  phenomena  happen.  Different  from  Case  I,  the  variation  of  L  and  Ksym  decrease  obviously.  It  means  that  L  and 
K sym  vary  a  little  when  the  parameters  are  changed  in  these  three  cases  and  they  are  the  smallest  in  non-linear  model. 
In  conclusion,  the  density  dependence  of  coupling  constants  can  make  L  and  Ksym  sensitive  to  the  changing  of 
parameters,  and  Case  I  is  most  likely  to  reproduce  the  experimental  results  for  it  being  the  most  sensitive. 


L/MeV 


L/MeV 


Figure  5:  The  Same  as  Figure  3  at  0.5, 1,  2  and  3  Times  Saturation  Density  for  Comparison 
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4.  CONCLUSIONS 

In  the  mean  field  approximation  of  nonlinear  relativistic  a-CO-p  model,  we  study  the  influence  of 
density-dependent  coupling  constants  on  binding  energy  of  a  nucleon  and  S(pB )  of  infinite  symmetric  nuclear  matter. 
We  found  that  greater  values  of  Tp  lead  to  greater  Ksym  and  L  when  Ta  c  and  d  are  fixed  and  indicate  stiffer 
S(pB)  in  high  density  region.  In  addition,  the  density  dependence  of  coupling  constants  can  not  only  decrease  the 
symmetry  energy  in  higher  density  region  but  also  make  L  and  Ksym  sensitive  to  the  changing  of  parameters  and  Case  I  is 
most  likely  to  reproduce  the  experimental  results  for  its  highest  sensitivity. 
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